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We develop an analytical approach for probability amplitudes of Kapitza-Dirac effect that merge 
together the Raman-Nath and Bragg regimes of interaction. 
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The Kapitza-Dirac effect constitutes diffraction of a structureless particle (electron) by 
a standing electromagnetic wave, exhibiting the particle-wave dual nature of matter in one 
of most convenient ways [l|. It appears as a counterpart to familiar optical diffraction 
by the periodic material grating. Matter wave of the particle beam plays the role of the 
incoming wave while the spatial periodicity of the grating interaction is ensured by the 
periodic structure of the optical standing wave potential. Since its prediction in 1933, this 
phenomenon was addressed theoretically many times [2|, y] and has b_een nicely observed 



experimentally by H. Batelaan et. al at Nebraska- 



a 



incoln University |^. (Detailed content 



can be found in review article J5| and dissertation 

In frame of ID model of sinusoidal periodic potential the electron wave function has 

a form 

+00 

■q/(z,t)= ^ Cn{t)exp[i{no + n)2kz] (1) 

and the problem reduces to the following difference-differential equation for c„ (t) probability 
amplitudes of the diffraction modes|7|: 



+ za;^(rao + n) + z?7o 1 c„ (t) = -i— - (c„_i (t) + c„+i (t)) , (2) 

n = 0, ±1 ± 2, .... Here k is the wave number of the running waves which constitute the 
periodic potential, no = Pinuiai/'^^k stands for the normalized electron initial momentum, 
ujr = Er/h with Er = {2hk) /2m is the so called recoil frequency and Uq represents the 
amplitude of the ponderomotive potential in h units. In representation (1) the momentum 
transfer occurs in discrete units 2hk and generally is interpreted as absorption/stimulated 
emission of photon pairs from the counterpropagating travelling waves, which give the stand- 
ing wave. 

To clearly identify the limiting Raman-Nath and Bragg regimes it is convenient to intro- 
duce a new amplitude 

Cn if) = z" exp [iUot + iur {no + n) t\ Cn it) , (3) 

which transforms Eq.(2) to 

^Cn it) = ^ (e'"'-(2-o+2-i)*C'„_i(t) - e-*"'-(2"o+2"+i)*C'„+i(t)) . (4) 



The analytic solution, presenting the Raman-Nath regime of interaction, corresponds to the 
limiting case Ur {2no + 2n — 1) t << 27r, when the system (4) loses the time- dependent 
exponential coefficients, transforming into equation with ffist kind Bessel function solution: 
Cn{t) = Jn{U{)t). Bessel function population flow of diffraction modes has a dominantly 
double-peaked pattern, symmetrically distributed about the initial state . This regime of 
realization assumes that the change of the electron position along the standing wave direction 
changes negligibly compared to the standing wave spatial period. 

The second, Bragg regime of interaction distinguishes only discrete, equidistant values 
or the electron initial momentum, namely in our notations n^ = ±1/2 ± 1, ±3/2.... Taking, 
for example, condition for the first order diffraction tiq = —1/2 {po = —hk), one can easily 
see that the two amplitudes at the right hand side of system (4), Co{t) and Ci(t), lose the 
time dependence in exponential coefficients, while the other ones preserve it. Assuming now 
an additional condition that Ur (2no ± 2n — 1) t >> 27r for any Ur and n (except, of course, 
n = 0, 1), we get rapidly oscillating coefficients for terms n ^ 0,1 and thus almost totally 
suppress their contribution to the final result (it reminds the rotating wave approximation 
widely used in the theory of matter-laser resonance interactions). After neglecting all these 
terms, one arrives to a simple pair of equations 

|co(t) = -^Ci(t), |ci(t) = ^Co(t), (5) 

resulting in Co(t) = cos{Uot/2) and Ci(t) = sm{Uot/2) probability amplitudes for direct 
{n = 0) and Bragg diffraction (n = 1). 

There is no exact analytical solution to Kapitza-Dirac problem in frame of Schrodinger 
equation, which will be valid for any interaction time periods and free of strict limitations on 
the system parameters. Diffraction regularities have analytically been treated in mentioned 
regimes of interaction and in close neighborhoods. They favor short- and long-time regimes 
respectively and are also known as the thin- and thick-crystal approximations. 

In this paper we develop a theory which treats the quantum particle (electron) diffraction 
in the ID periodic potential for any times of interaction. Our formula quantitatively correctly 
describes both Raman-Nath and Bragg regimes of interaction, thereby merging them into 
one essence of diffraction process. 

Our approach to the infinite system of equations (2) originates from the remark that it 
connects the seeking amplitudes with opposite parity on the left hand and right hand sides 



of the equations. Initial condition c„(0) = Sn,o iSj however, different for these two famihes 
of amphtudes: all odd ones are zeroes, while the only nonzero member sits in the even- 
manifold. Any approximate solution should be sensitive to the initial conditions too, as 
we get rid of one of the parities in the set of equations by introducing new, phase-shifted 
amplitudes 

c;r(t)=i"exp[if/ot]c„(t) (6) 

and arriving to 

(I + ioor ((no + n)^ + 1))^ (I + iujr (no + nf) c;;{t) + 
Auol (no + nf (| + loOr (no + nf) + f (| + loOr ((no + nf + l))) c;i{t) 



^ I -IT +^Wr ((no + n) + Ij I (C„_2 



~{t) + c;;^{t)). 



It preserves the original tridiogonal form of (2) but is now a third order differential - difference 
equation. In the following treatment the set (7) will be regarded to describe the even- 
manifold of amplitudes. 

We look for a trial solution of a definite integral form 



.^\n r^ 



Cn(t) = exp [—iur (no + n)^ t] / cos (nip) exp [iXn{(p)t] dip, (8) 



TT JO 

where the n-dependent function A„((/9) has to be determined yet. Inserting Eq.(8r into (7) 
we obtain the following three difference algebraic equations: 

An-2(v^) +4a;r (no + n- 1) - A„(v9) = (9) 

A„+2(v^)+4wr(no + n+l)-A„(v9) = (10) 

(A„(v9) + Urf A„(v9) - 4w^ (no + nf A„(v9) - -^ (An(v^) + ^r) j COS (rup) (11) 

f/2 f/2 

= -J- (A„-2(v') + Wr (6no + 6n - 3)) cos ((n - 2)p) + -j- (A„+2(v^) - ^r (2no + 2n - 3)) cos ((n + 2)(p) 

Hence,for the trial function to be an exact solution of Eq.(7), the Eqs. (9)-(ll) should be 
identical to each other. Eq. (9) and (10) are really mutually equivalent and one of them, 
for instance Eq. (9), can be put out from consideration. The last one, Eq. (11), however, 
is not equivalent to Eqs. (9) and (10). This means that the analytic form (8) can not be 
an exact solution to the problem. Our approximation lies just in this point and amounts 



to taking of Eqs. (9)-(ll) as a system of three equations relative to A„+2(v^), ^n-2{f) and 
seeking A„ (</?). Then finding A„+2(v^) as a hnear function of A„(v?) and inserting it into Eq. 
(11), we arrive to a third order algebraic equation 

(A„(y^) + 2ur/3f + pn if) (Xnif) + 2u;./3) + g„ {^) = (12) 

with real coefficients 

Pn iv) = - (f/o cos^ {lp) + Aul {no + nf + u^/S) 

and 

Qn (V^) = -^^r - UoUJr f - COS^ {^) + (?^0 + "-) COS {2ip) j + -U^ (% + uf . 

The sign of the first coefficient is always negative and depending on the sign of quantity 



.37 V 2 
one has two distinct forms for the solutions of Eq. (12) [8J. 

Later on we'll denote the three roots of Eq. (12) as A„(j; <f), j = 1, 2, 3, the corresponding 
amplitudes as c^{j; t) and the respective wave functions as \E' (j; z, t). Then the general form 
of probability amplitudes should be written as 

^ {z, t) = /ii^ (1; z, t) + h2^ (2; z, t) + h^^ (3; z, t) (13) 

or, equivalently 

( — 1)" 
Cn(t) = exp \-iUr iuQ + nft- iUqI] ^x (14) 

L J ^ 

/ {hi exp [iA„(l; (p)t\ + /i2 exp [iA„(2; (p)t\ + /13 exp [iA„(3; (■p)t]) cos {rup) dip. 
Jo 

Three /i-coefficients are determined from initial conditions for the wave function (13) and 

its first and second derivatives: 

hi + h2 + hs = 1, 

hh + l2h2 + hh = 0, (15) 

Lihi + L2h2 + Lsh3 = U^/2, 

with the following notations for coefficients: 

^j = - K=o{J;'p)d!p, Lj = - Xn=o{j;<ffd!p, (16) 
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FIG. 1: (Color online) n = diffraction probability amplitude evalution in Bragg diffraction 
regime. Solid line gives the result of exact numerical simulation. Dashed line is the amplitude 
evolution for analytic formula described in the text. The dashed line is shifted vertically by 
0.1 to make the two lines visually different (we will make this shift in all figures) . Horizaontal 
time axis is normalized to green light {Xgreen = 6000A) Compton backscattering frequency shift 
^^Comptori = '^'^fi^green/ ^eC^ , where Mg is the cMctron mass. Other parameters are no = —0.5, 

Uq = 0. 01 Auj Compton/'^ , '^r = ^^Comptonh- 

j = 1, 2, 3. This step crowns the procedure and hence the formula (14) presents the solution 
of the problem for any even n, the amount of acquired momentum in 2fik units . 

To determine still untouched odd- probability amplitudes, we have to return to original 
equation (1) and shift the numbering by one: 

— + iujr {no + nf + iUo j c„+i (t) = -i^ (c„ (t) + c„+2 {t)) , (17) 

Inserting even-ra solutions into the right hand side of (17) and simply integrating equation 
with zero initial condition, we complete our approach to analytic solution of the stated 
Kapitza-Dirac diffraction problem. 

To value the developed analytic approximation, we have compared its results with the 
exact numerical solutions of the original set of Eqs. (2). In order to implement these simu- 
lations we have used the Crank-Nicolson method [9|. Comparison shows that the presented 
approximation works excellent in both, Raman-Nath and Bragg regimes of interaction. Two 
particular cases are illustrated in Figs. 1-4. The graphs in each figure are indistinguishable 
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FIG. 2: (Color online) n = 1 diffraction probability amplitude evalution in Bragg diffraction regime. 
All the parameters are as in Fig.l. 
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FIG. 3: (Color online) n = diffraction probability amplitude evalution in Raman-Nath diffraction 
regime, no = -0.5, Uq = Aujcomptonf-n- , ujr = 0. 001 Aujcompton/-^- 



from each other at sight and are shifted in vertical direction in illustrative purposes. 

In intermediate regimes (relative to Raman-Nath and Bragg) our numerical calculation 
has definite restrictions, connected with the limitations on calculation of inverse matrices 
required by the Crank-Nicolson method. Thus ultimate conclusions here cant be done 
yet. However, every case when we were sure in the correctness of numerical calculations 
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FIG. 4: (Color online) n = 2 diffraction probability amplitude evalution in Ranian-Nath diffraction 
regime. All the parameters are as in Fig. 3. 
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FIG. 5: (Color online) n = diffraction probability amplitude evalution in intermediate diffraction 

regime, no = -0.5, Uq = Aucompton/^T^ , ^r = AoJCompton/T^- 



the coincidence between our approximate analytical results and numerical ones was quiet 
good. Fig. 5 and 6 illustrate such a case with parameter values Uq ~ Ur (note that the 
Raman-Nath approximation requires Uq » Ur, and the Bragg approximation - the opposite 
one Uq « Ur ). 

This gives some credibility to the presented analytical approximation over the intermedi- 
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FIG. 6: (Color online) n = 1 diffraction probability amplitude evalution in intermediate diffraction 
regime, no = -0.5, Uq = Aujcompton/'i'n- , uJr = Aujcompton/T^- 

ate range of parameters too and in the future we will endeavour in reaching a full definiteness 
in this direction too. 

This work is supported by Alexander von Humboldt foundation and NFSAT/CRDF 
Grant No. UCEP-0702. 
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